UNIFORM BIRKHOFF INTERPOLATION WITH RECTANGULAR 

SETS OF NODES 



MARIUS CRAINIC AND NICOLAE CRAINIC 

Abstract. In this paper we initiate the study of Birkhoff interpolation schemes with 
main emphasis on the shape of the set of nodes. We concentrate here on the simplest 
shapes ("rectangular"). The ultimate goal is to obtain a geometrical understanding of 
the solvability. We present several regularity criteria, a few conjectures, with geometric 
interpretations, and we illustrate with many examples. 
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1. Introduction 

Apart from its historical importance (see e.g. Newton's interpolation formula, Gauss' 
quadrature formula), interpolation by polynomials still plays a central role in the local 
construction of various approximation schemes (splines, finite elements, cubature formulas, 
etc). The homogeneous Birkhoff- Hermite problem is one of the most general multivariate 
polynomial interpolation problems. In general, it depends on a set Z C W 1 of nodes, a 
"lower set" S defining the "interpolation space" Vs of polynomials (to which the solution 
is required to belong to), and a set A C Z™ of derivatives which appear in the interpolation 
equations. The univariate case (n = 1) is quite well understood, and it behaves fundamen- 
tally different from the multivariate case (Hj. For notational simplicity, we restrict here the 
the bivariate case. 

In the literature one finds two types of results. On one hand, there are criteria whose 
conclusion holds in the generic case (when the points of Z are in general position). A 
very good example is the work of R.A. Lorentz On the other hand, there are more 
constructive results which use ideas from the univariate case to produce explicit solutions 
for some special classes of interpolation schemes. A good example is the work of Gasca and 
Maeztu [7] (for more complete list of references, and historical comments, we refer to the 
same |Z1|!5|). I n such results, although the shape of Z is very special, it is very rare that it 
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is the starting point. More precisely, the shape comes on the second place: it is the one for 
which the method works. 

Clearly, the patterns of the general problem are far too complex to be understood in an 
unified manner. One of the few common features of all approaches is that the shape of Z 
(and sometimes that of A too) plays an essential role. On the other hand, when it comes to 
applications, it often happens that the set of nodes Z is given, has a very particular shape 
(think e.g. of cubature formulas), and there are only special cases/aspects of the problem 
that need to be solved/understood. 

In this paper we initiate the study of multivariate Birkhoff interpolation schemes with 
main emphasis on the shape of Z. After a short discussion on "shapes" , we will concentrate 
on the simplest shapes, namely the rectangular ones fsee 12.5(1 : 

Z = {(x i ,y j ):0<i<p,0<j<q} 

(p and q are non-negative integers). Apart from the main definitions, and a short discussion 
on "cartesian shapes" and their relation with the uniqueness of Birkhoff-Lagrange schemes 
(see !2.4|) . section El also brings together the main constructions that are relevant to the study 
of rectangular sets of nodes (most importantly, the notion of blow-up presented in I2.3[) . At 
the end of the section we briefly compare the generic case with the rectangular case. In 
section we present several regularity criteria which are very useful in examples, and we 
derive several consequences. Inspired by all the results and all the examples, we conjecture 
that lower sets S that are part of a regular scheme with rectangular sets of nodes must 
be the result of a blow-up, and we relate this to a stronger version (for rectangular sets of 
nodes) of the Polya inequalities. The results presented in this section can also be viewed as 
a confirmation of the conjectures in many unrelated cases (depending on the size and shape 
of A). At the end, we prove that the conjectures hold true in one more case: p = q = 1 
(and, this time, for all A). On the other hand, by looking at the pictures associated to 
regular schemes, one observes another stricking property that regularity seems to imply (a 
geometric property this time!). Initially, it was tempting to formulate this property into 
yet another conjecture, but, to our surpprise, it turned out to be equivalent to the the 
Polya-type inequalities mentioned above. This is the content of Theorem 14.11 in Section |IJ 
In the same section we point out yet another geometric property that regular schemes seem 
to share (Conjecture 01 which we use as a guide for constructing interesting examples, cf. 
e.g. Example I5.10JI . And, finally, Section contains a large list of examples. 

Finally, we would like to point out that, although the results presented here are quite 
satisfactory for handeling large classes of examples, we feel that the main achievement 
of this paper is the understanding we gain (e.g. that blow-ups are relevant, that A is 
geometrically related to S, etc). In this direction, we also believe that finding relations with 
other fields (e.g. algebraic geometry or algebraic topology) would give a new understanding 
of the conjectures above and of the interpolation problem (e.g., note that "the degrees of 
multivariate polynomials" together form an operad, and the blow-up is one of the simplest 
operations associated to it; see 12 .3(1 . The existence of a deeper such relation is (yet) another 
conjecture we would like to adress here. 

This paper is part of the second author's PhD dissertation, and it has circulated as 
a Utrecht University preprint. The present version follows "qualified suggestions" and is 
considerably shorter (in particular, it skips some of the technical but completely elementary 
proofs which will be presented elsewhere jH HI El HI IS] ) • 
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2. Main definitions and constructions 

2.1. Main definitions/terminology: In this paper we will use the following notations/terminology: 

• Z denotes a finite set of points in M 2 , which plays the role of the set of nodes of the 
interpolation problem. 

• A denotes a finite set of pairs of non-negative integers, describing the (order of 
the) derivatives appearing in the interpolation problem. Typical examples are the 
rectangles R(u,v) and the triangles T(n): 

R(u,v) = eZ 2 + :i<u,j< v}, 

T(n) = {(i,j) eZ 2 + :i+j <n}. 

• S is a lower set, that is, a finite subset Sc^ with the property that 

(u, v ) <G S => R(u, v) C S. 
Corresponding to S is the space 

V s = {PeR[x,y]:P= £ a M xV}, 

(i,j)es 

and Vs plays the role of the interpolation space (where the interpolation is taking 
place). 

• such triples (Z, A, S) are called interpolation schemes, and the associated problem 
is: for given constants {cij{z) : (i,j) € A, z G Z}, find polynomials P £ Vs 
satisfying the equations 

Qi+jp 

(i) ^W (z) = c ^' yzeZ > ( ^' )GA 

One says that (Z, A, S) is solvable if the problem has solutions for any choice of the 
constants. One says that (Z, A, S) is regular if the problem has unique solution. 
The simplest necessary condition for regularity is normality of (Z, A, S), that is, 
l^l = All schemes in this paper are assumed to be normal. 

• the interpolation equations are in fact a system of linear equations, and one denotes 
by D(Z, A, S) its determinant (well defined because of the normality condition). 
This is a polynomial on the coordinates of the nodes (2n such coordinates, where 
n = \Z\). Hence, fixing (A,S), if (Z, A, S) is regular for one choice of Z, then it 
is regular for almost all choices of Z (where "almost all" refers to the Lebesgue 
measure on M 2n ). In this case one says that (A, S) is almost regular with respect to 
sets of n nodes. 

2.2. Representing lower sets: First of all, for any non-negative integers s, no > n\ > 
...>n s , 

S y (n , . . . ,n s ) = {(i,j) £Z 2 :fl<i<s,0<j< n{\ 

is clearly a lower set. Conversely, any lower set S can be uniquely written in such a way. 
Given any set B C Z 2 , we introduce the notation 

B y [a] = {0: (a,0) eB}. 

Then, given a lower set S, the associated s is the maximal number i with the property that 
S y [i] ^ 0, while ni is determined by <Sy[i]: 

S y \i] = {0,l,...,ni} 
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A similar discussion is obtained by interchanging the role of x and y. Representing S both 
as Sy(nQ, . . . , n s ) and as S x (mo, mi, . . . , m t ), the relation between the two is: 

t = no, nij = max{i : rii > j}. 




S y (9, 9, 9, 7, 7, 7, 5, 2, 2, 2, 2)= ^(10, 10, 10, 6, 6, 6, 5, 5, 2, 2) 



Figure 1. 

Another way of representing lower sets is by making use of its "exterior boundary points" . 
Given S, a point (u, v) <G S is called a boundary point if (u + 1, w + 1) ^ S\ We denote by 
35 the set of such points. We consider the following two possibilities: 

(i) (u,v + l)€S; 

(ii) (u+l,v) G 5. 




exterior boundary interior boundary x-direction y direction 

point point boundary point boundary point 



Figure 2. Boundary points. 

We denote by d e S the set of boundary points (u, v) for which neither of the two conditions 
above is satisfied ("exterior boundary points"), by diS the set of those which satisfy both 
conditions ("interior boundary points"), by d x S the set of those for which only (ii) holds 
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true ( "x-direction boundary points"), and by d y S the set of those for which only (i) holds 
true ( "y-direction boundary points"). These four sets form a partition of the boundary 
d(S) of S. 

The set d e S of exterior boundary points determines S uniquely, since 

S= [J R(u,v). 

(u,v)ed e s 

It is actually useful to label the elements of d e S by (a^, 61), (cifc-i, 62), . . . , (a±,bk), indexed 
so that a\ < a>2 < ■ ■ ■ < Ofc, and b\ < 62 < . . . < b k . 




Figure 3. Exterior boundary points. 

Then an element (u, v) is in S if and only if 

ai-i <u<ai, <v < bk-i+i 

for some / € {0, 1, . . . k}. Also, 

= {(a,, bj) :i + j = k + 1}, diL = {(en, bj) :i + j = k}, 
d x S = {(u, b k _j) : 1 < j < k,aj < u < a j+1 }, 
d y S = {(a k -i, v) : 1 <i < k,bi < v < b i+1 }. 

2.3. Blowing up lowers sets: Central to Birkhoff interpolation with rectangular sets 
of nodes is the notion of blow-up. Given p, q non-negative integers, and a lower set S, we 
definer a new lower set S p ' q which is obtained by "blowing up" each of the points of S to 
a (p, c/)-rectangle. One may think that a copy of R(p, q) is placed on each of the points 
of S, and then one pushes these rectangles minimally to the right and upwards until they 
become disjoint. For an example with p = 1, q = 2, see the picture. 
More formally, 

SP.« = {(a,/3):([^ I ],[^ I ])6 5}. 
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Figure 4. Blowing up. 

If S = S y (riQ, . . . , n s ), then S p ' q = S y (n , . . . , n' s ,), where 

k' = {q + l)(k + 1) - 1, n'i = (p + l)(nj + 1) - 1) where j = [— ^— 1. 

q + 1 

In terms of the exterior boundary points, d e S p ' q consists of the pairs ((p + l)u, (q + l)v) 
with (u, v) £ d e S. 

It is interesting to point out the structure underlying these operations, structure that 
becomes even more important in higher dimensions. Let us denote by £{n) the set of all 
lower sets L C . Then, for any positive integers n, k%, . . . , k n , there are "operations" 

C(n) x C(ki) x ... x C(k n ) — > C(ki + . . . + k n ) 

as follows: given (S,Si, . . . , S n ), the result S' = S • (Si, . . . , S n ) of this operation is the 
lower set with the property that 

T s > = {P{Pi(xi, . . .,x kl ),P 2 (x kl+1 , ...,x kl+k2 ),...):PeVs,Pi£V Sz }. 

These operations make {C(n)} into an operad [S], which we will call the multi-degree op- 
erad. The reason for the terminology is that lower sets can be viewed as an extension to 
higher dimension of the ordered set of natural numbers (non- negative integers): any natural 
number n defines a lower set [n] = {0, . . . ,n}, and these are all the possible lower sets in 
dimension one. In general, an n-multivariate polynomials will have an associated lower set 
playing the role of degree, and Vs will be the space of polynomials "of degree (at most) 5"' . 
In particular, lower sets arise naturally also when discussing symbols of multi-deifferantial 
operators. Note also that the blow-up is one of the simplest operations encoded in the 
operad {£(n)}: 

s^ = s-(\p],[q})- 

Shapes more complicated then the rectangular one would probably require the use of the 
other operations associated to {C(n)}. 
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2.4. Particular shapes: We now give a more precise meaning to the "shape" of a set of 
nodes. First of all, given a lower set S, one says that a set of nodes Z is S-cartesian if it is 
of type 

Z = {(xi,Vj) : 0', j) G S}, 

where the x%s are distinct real numbers, and similarly the y/s. We say that Z is cartesian 
if it is S-cartesian for some lower set S. 

On the other hand, any set of nodes Z induces two lower sets S X (Z), and S y (Z), which 
reflect the shape of Z. To describe S y (Z), one covers Z by lines Iq, . . . , Ik parallel to the 
OY axis, and define the numbers so that on each line l( there are exactly m + 1 points of 
Z. We index the lines so that no > n\ > . . . > n^, and we define S y (Z) = S y (no, . . . , n/t). 
The lower set S X (Z) is defined similarly, by interchanging the role of x and y. 

The relation between the two is that a set of nodes Z is cartesian if and only if S X (Z) = 
S y (Z). In turn, this implies that cartesian sets of nodes show up naturally in relation 
with the uniqueness of Birkhoff-Lagrange schemes. Birkhoff-Lagrange schemes are those 
for which A contains only the origin (and is then suppressed from the notation). Given 
any Z, there exists at least one lower set S so that the Birkhoff-Lagrange scheme (Z, S) is 
regular. However, S is unique if and only if Z is cartesian. For details, see pQ. 

2.5. Rectangular shapes: Of special interest for us are the cartesian sets with respect to 
the rectangles R(p,q). Such sets are called (p, q) rectangular (or just rectangular if we do 
not want to emphasize p and q). Hence, Z is (p, q) rectangular precisely when it is obtained 
by intersecting (p + 1) vertical lines with (q + 1) horizontal lines. Writing 

z = {(x i ,y j ):0<i<p,0<j<q} 

where the distinct real numbers, and similarly yo,...,y q , the determinant 

D(Z, A, S) is a polynomial on only t = p + q + 2 variables (namely the ccj's, and the y/s). 
We then say that a pair (A, S) is almost regular with respect to (p, g)-rectangular sets of 
nodes if there exists such a Z so that (Z, A, S) is regular. As in the generic case, this implies 
that (Z, A, S) is regular for almost all choices of the (p, q) -rectangular set s Z (this time, of 
course, the "almost" part refers to the Lebesgue measure on M*). The regularity is defined 
similarly. Of course, one can also allow complex nodes, and talk about (almost) regularity 
with respect to (j>, q) -rectangular sets of complex nodes. 

2.6. Examples: To get a feeling about the effect that the shape of Z has on regularity, 
let us point out two cases. For A = {(0, 0), (1, 1), there are no regular schemes whith 
rectangular sets of nodes, while there are many regular schemes with the set of nodes in 
general position. Consider now the case where 

A = {(0,0), (0,1)}. 

Then one can show that all schemes (^4, S) in which S has the property that |S*| = 2n and 
S contains at most n elements on the OX axis, are almost regular with respect to sets of n 
nodes. On the other hand, there is only one scheme which is almost regular with respect to 
(p, q) -rectangular sets of nodes, namely S = R(l,0) p ' q = R(2p + l,q). These will become 
clear in the next sections. 

2.7. Polya conditions: To compare the multivariate case with the univariate one, and 
general sets of nodes with rectangular ones, it is interesting to look at Polya conditions. 
These are certain algebraic inequalities that are forced by regularity. The classical Polya 
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condition [5] say that, given the regular scheme (Z, A, S) then, for any lower subset IcS, 

(2) n\LnA\>\L\, 

where n = \Z\. This follows quite easily by counting the numbers of zero in the determinant 
associated to the interpolation equations Moreover, in the limit case (when equality 
holds) (Z,A n L,L) must be regular too. This corollary applies to the univariate case as 
well. Writing A = {ao, ax, ■ ■ ■ ,a s } with clq < a\ < . . . < a s , the Polya conditions become: 

di <n ■ i, V0<i<s. 

Moreover, this condition actually insures regularity for almost all sets of nodes Z 6 . More 
precisely, given (A, S) with \S\ = n\A\, (A,S) is almost regular if and only if it satisfies 
the Polya conditions. Moreover, if n = 2, then the Polya conditions are sufficient also for 
regularity. Such properties do not hold in the multivariate case. 

When the shape of Z is particular, one can strengthen these inequalities. For instance, 
using boundary points (see 12.21 above), and by a careful analysis of the zeros in the deter- 
minant associated to interpolation equations, one can prove [2] that if (Z, A, S) is regular 
and Z is (p, g)-rectangular, then, for any lower subset L C S, 

n\AC\L\ > \L\ +pq\ADdL\ + (p + q)\A n d e L\ +p\And y L\ + q\And x L\. 

3. Regularity criteria 

In this section we bring together several regularity criteria for uniform Birkhoff inter- 
polation schemes with rectangular sets of nodes, then we present some consequences and 
then point out a property that schemes with rectangular sets of nodes seem to satisfy (as 
suggested by the general criteria, proven for \A\ < 6, and also for p = q = 1). Since the 
proofs of these criteria, although a bit technical, are completely elementary, they will be 
presented elsewhere IU E] ■ Many examples illustrating the usefulness of the criteria will 
be presented in the last section. 

Throughout this section, (Z, A, S) is a uniform Birkhoff scheme, with a (p, </)-rectangular 
set of nodes 

Z = {(x l ,y J ):0<i<p,0<j<q}. 
We denote by A x the elements of A on the OX axis, we also use the similar notations A y , 
S x , S y , and we put Z x = {xq, . . . , x p } (the projection of Z on the OX axis), and similarly 

Zy. 

Proposition 3.1. If (Z, A, S) is regular, then 

\A\ < | A x | I Ay | . 

Moreover, if the equality holds, then S must be: 

S = R(p',q'), p' = (p + l)\A x \-l, ^ = (9 + 1)1^,1-1. 

This property puts restrictions on the set A. For instance, if mix(A) is the number of 
mixed derivatives coming from A, one immediately gets 

(3) y/mix(A) < V\A\-1. 

The next property takes care of the case where one of the lower sets 5*^(^4) or S X (A) 
associated to A (see I2.4j) is rectangular. We then have (for the notation A, [a], see 12 .2|) : 

Theorem 3.2. If Z A has the property that S y {A) is (s,t)- rectangular (in particular, if 
A is (s,t)- rectangular), then the scheme (Z,A,S) is regular if and only if 
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(i) S = R(p',q'), with p' = (s + l)(p + 1) -l,q' = (t+ l)(q + 1) - 1. 

(ii) The univariate scheme [Z x , A x , S x ) is regular. 

(iii) All the univariate schemes (Z y ,A y [a],S y ), with a € A x , are regular. 

The next criterion shows the influence that lower subsets of A have on S. 

Theorem 3.3. If (Z, A, S) is solvable (i.e. the interpolation equations have at least one 
solution, but not necessarily unique), and if A contains a lower set R, then S must contain 
the lower set R p ' q . 

In particular, if A is lower, one can completely clarify the situation. The outcome makes 
use of the univariate Hermite polynomials for the set of nodes Z x . For non-negative integers 
u, a, and a node x s , we consider 

ul(k — u)\ m s 

k=u 

where 

<i> a {x) = (x- x Q r +i ...(x . . . (x - Xp r +i . 

We make the convention that s (x) = if a < u. Similarly, we denote by H?'* the 
univariate Hermite polynomials associated to the set of nodes Z y . We also write 

d e S = {(a k ,h), (ai,b k )}, 

as in 12.21 With these, we have: 

Theorem 3.4. If the set A of derivatives is lower, then (Z, A, S) is regular if and only if 
S = A^ q . 

Moreover, in this case, for any (u,v) € A, and any node (x s ,yt) G Z, the polynomial 

Huf v (x,y) = 

H^ s (x)(H b v *\y) - H h v ^\y)) + H^ s (x)(H b v k ^\y) - H b v k -' 2 \y)) + 

. . . + H a u k -^ s {x){H b v ^{y) - H b ^(y)) + H?* \x)H^' x \y) 

is the fundamental interpolation polynomial at the node (xi,yj), with respect to the deriva- 
tive (u, v), i.e. 

d a+l3 H^ v n f 1 ifa = u,P = v,i = s,j = t 
dx a dyP ^ Xl > y i) ~ \ otherwise 

The next property determines S "around the coordinate axes" . We use again the nota- 
tions Sy[f3], S x [a] introduced in 12.21 

Theorem 3.5. If (Z, A, S) is a regular UR Birkhoff scheme, where the set of nodes Z is 
(PiQ) -rectangular, then 

S x [0) = S x [l] = ... = S x [q) = {0,1,..., p'}, 

Sy[0]=S y [l] = ... = S y [p] = {0,1,..., q'}, 

where p' = {p+l)\A x \ - 1, q' = (q + l)\A y \ - 1. 

And, finally, the following shows how one can move or remove points of A on the coordi- 
nate axes. We denote by S' the set obtained from S by removing the lowest {p, g)-rectangle 
on its most right (that is, remove the last (p + 1) elements of S from each of the lines y = i, 
with < i < q). Note that S' is lower if and only if \S x (q + 1)| < p' — (p + 1). 
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Theorem 3.6. Let (Z, A, S) be an uniform Birkhojf scheme with rectangular set of nodes. 
Then: 

(i) If (Z, A, S) is regular, then the univariate scheme (Z x , A x , S x ) must be regular. 

(ii) If A C Z 2 is obtained from A by moving some of its elements on the OX axis, then 
(Z, A, S) is still regular if (Z x , A x , S x ) is. 

(iii) If A' is obtained from A by removing one of its elements from the OX axis, S' is 
lower, and and both univariate schemes (Z x , A x , S x ) and (Z x , A' x , S' x ) are regular, 
then (Z, A, S) is regular if and only if (Z, A' , S') is. 

There is an obvious variation obtained by interchanging x and y. Also, there is a par- 
ticularly good way of changing A. Let us denoted by A max the set obtained by replacing 
the elements of A on the OX axis with {0,p + 1, 2(p + 1), . . . , (p + 1){\A X \ - 1)}. This is 
the "maximal replacement" allowed by the univariate Polya condition fsee l2.7|) . Moreover, 
(A ma x,x, S x ) is automatically regular with respect to all sets of (p + 1) nodes (under the 
normality condition \S X \ = (p+ l^vl^l). In particular, applying first (ii) and then (iii), we 
deduce that (Z, A, S) is regular if and only if (Z x , A x , S x ) and (Z, A' max , S') are. 

A large number of examples which show the usefulness of these criteria are presented 
in Section |31 Let us now point out several consequences. First of all, using Theorem 13.51 
to determine S, and then Theorem I3.nl to reduce the size of the problem (and A max to 
simplify), we immediately deduce 

Corollary 3.7. If A contains no mixed derivative, then (Z, A, S) is regular if and only if 
the univariate schemes (Z x , A x , S x ) and (Z y , A y , S y ) are regular and 

S = T p '", where T = R{\A X \ - 1, 0) U R(0, \A y \ - 1). 

The next case (one mixed derivative) is easier to state in terms of almost regularity 
instead of regularity. 

Corollary 3.8. Assume that A contains only one mixed derivative, call it (a, (5), and let 
{(aj,0) : < i < s} be the elements of A on the OX axis, and {(0,bj) : < j < t} the 
ones on the OY axis. Then (A,S) is almost regular with respect to (p,q) -rectangular sets 
of nodes if and only if s,t > 1, 

<H <»(p + l),&j <j(q + l),a<2p + l,P<2q + l 

(for all i and j), and 

S = T p ' q , where T = R(\A X \ — 1, 0) U R(0, \ A y \ — 1) U 1). 

Proof. The conditions on the a^'s and the bj's are just the univariate Polya conditions 
(see 12. 7(1 . On the other hand, if we require that no nontrivial polynomial of type P = 
(y ~ 2/0 ) • • • {y ~ yq)Q{ x ) ^ 'Ps ma y satisfy the homogeneous interpolation equations of the 
scheme, we see that \S x [q + 1]| < 2p + 1, and, similarly \S x \p + 1]| < 2q + 1. Using also 
Theorem 13. 51 and the fact that S is lower, these inequalities immediately imply the last part 
of the statement. This allows us to remove elements of A from the coordinate axes (use 
Theorem 13.61 and A max ) and reduce the problem to s = 1 and t = l,a\ = (p+1), b% = (q+1). 
Finally, if the inequalities on a and (3 do not hold, we have (a,/3) 3 R(p + l,q+ 1), and, 
using the Polya inequality for L = R(p + 1, q + 1) (cf. 12.7(1 we would obtain 

3n > (n + p + q + 3) + 2pq +p + q, 

i.e. 3n > 3n + 1 which is impossible. Conversely, if a and satisfy these inequalities, then 
we can move (p+1, 0) G A to (a, 0), and then S y (A) becomes rectangular and we can apply 
Theorem □ 
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Corollary 3.9. If(A,S) is almost regular with respect to (p,q) -rectangular sets of nodes, 
and \A\ < 6, then S is of type T p,q for some lower set T . 

Proof. Using ©, A will contain at most one mixed derivative except for the case where 
\A\ = 6 and A contains two mixed derivatives. The first cases follow from the previous 
corollary, while the last one follows from the limit case of Proposition 13.11 □ 

The results above (see also below, and the section on examples) suggest the following 

Conjecture 1. If a scheme (A,S) is almost regular with respect to (p,q) -rectangular sets 
of nodes, then S = R p,q for some lower set R. 

Let us give a stronger version of this conjecture, which also has a geometric interpretation 
that is easier to visualize (see the next section). Given a lower set S, we denote by n P:q (S) 
the number of elements (a,/3) £ S with the property that a is divisible by (p + 1), and 
f3 is divisible by (q + 1). The following can be interpreted as a Polya-type inequality for 
rectangular sets of nodes. 

Conjecture 2. (p,q)-Polya conjecture) If a scheme (A,S) is almost regular with respect 
to (p,q) -rectangular sets of nodes, then 

(4) \LHA\ >n m (L) 

for all lower sets L C S. 

Let us point out the relation with the Polya condition, as well as the relation between 
the two conjectures. 

Proposition 3.10. Consider the scheme (A,S) with \S\ = n\A\, where n = (p + l)(q + l). 

(i) For all lower sets L, n P:q (L) > —\L\. In particular, Conjecture^ is a strengthening 
of the Polya conditions HJ). 

(ii) S = R p,q for some lower set R if and only if ' n p ^ q {S) < \A\ (and then equality must 
hold). In particular, Conjecture^ implies Conjecture^ 

Proof. Given a lower set L, we define a new lower set L Ptq obtained from L by "collapsing 
(p, g)-rectangles to points" : 

L p , q = {(a, 0) : ((p + l)a, (q + l)/3) G L}. 

First of all, it is clear that n Ptq (S) = \S p>q \. Secondly, it is easy to see that (L p,q ) p q = L, 
and L C (L P}q ) p ' q for all lower sets L. Passing to cardinalities in the last inequality, we 
obtain (i). 

On the other hand, one immediately sees that S Ptq is the smallest lower set with the 
property its (p, g)-blow up contains S. In particular, if S = R p ' q for some lower set R, 
then R must coincide with S pq . So, the condition that S = RP' q for some lower set R is 
equivalent to saying that 5 C {S p ^ q ) p ' q becomes equality. Passing to cardinalities, we have 
\A\ < ripq(S), and the requirement is that equality holds. This proves (ii). □ 

Regarding (the stronger) Conjecture 12 it is verified in all examples we looked at, and 
the results of this section can be used to verify it in many other cases. Let us point out 
that it is true also when p < 1, q < 1. We state here the case p = q = 1, which will be very 
useful when looking at examples. 

Theorem 3.11. Given A and S , p = q = 1, the following are equivalent: 

(i) The scheme (A, S) is regular for some (1, 1) -rectangular set of nodes Z. 
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(ii) The scheme (Z,A,S) is regular for any (1,1) -rectangular set of nodes Z. 
In this case moreover, (A,S) satisfies Conjecture^ (hence also Conjecture^). 

Proof. By changing coordinates 

(x,y) i-» (2 1,2 1, 

x\ -x Q yi- y 

we can always assume 

Z= {(-1,-1), (-1,1), (1,-1), (1,1)}, 

and this proves the equivalence of (i) and (ii). Assume now that there exists a lower subset 
L C S such that nx t \{L) > \A fl L\. This precisely mens that R = Li,i (where we use 
the notations from the previous proof) satisfies \R\ > \A D L\. We now consider the space 
Vir of polynomials spanned by x 2a y 2 ^ with (a, 0) G R. Note that Vir C Vs- Next, since 
the polynomials in V2R depend on \R\ > \A\ variables (the coefficients), we find a non-zero 
polynomial P G V2R such that 

W J^(1,1)=0, V (M)^ 

But P G V2R implies that P(—x,—y) = P(—x,y) = P(x,—y) = P(x,y). Taking deriva- 
tives, these relations still hold true, at least up to a sign. In particular, (JSJ) implies that P 
satisfies all the equations Q corresponding to derivatives coming from Af)L, and with the 
constants Cij(z) equal to zero (the homogeneous equations). On the other hand, since L is 
lower and P G Vl, all the derivatives of P coming from A \ L are identically zero. Hence 
P G V2R C Vs would be a non-trivial solution of the homogeneous equations associated to 
our scheme, which contradicts regularity. □ 

Note also that, with the same argument as above (using the (p + l) th and the (q + l) th 
roots of unity) proves the following 

Corollary 3.12. If (A, S) is regular with respect to (p, q) -rectangular sets of complex nodes, 
then it satisfies Conjecture^ (hence also of Conjecture^. 

4. Geometric aspects 

In this section we discuss the geometric interpretations of the inequalities Q, which are 
easy to visualize on the picture (see the next section). 

Roughly speaking, this property says that S is obtainable by attaching a copy of the 
rectangle R(p, q) to each point of A, and then moving these rectangles upwards or to the 
right until they become disjoint. To make this more precise, we need to use "shifts" of A in 
S, which are transformations which move A step by step, at each step one of its elements 
being moved upwards or to the right on a new position which is still in S and which is not 
occupied by any other element of A. We also consider the set 

Z p ' q = {({p + l)i, (q + l)j) : i,j — positive integers }, 

which can be viewed as the blow up of the lattice Z of integral points situated in the first 
quadrant. 

Theorem 4.1. Given (A, S) (\S\ = n\A\, n = (jp+ l)(q + 1) ), the following are equivalent 

(i) \L n A\ > n Pt g(L) for all lower sets L C S. 

(ii) there exists a shift of A in S which moves A to S n Z p ' q . 
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Proof. Before starting the proof, let us fix some notations. Given an element e G S, 
we denote by e r the one step translation of e to the right, and similarly e u (upwards 

o 

translation). Also, for any lower set L, we consider two new lower sets L and L such that 

o o 

LC L C L, n Pjq (L) = n PA (L) = n p>q (L), 

o 

and such that L is the smallest possible one, while L is the largest possible. Note that 
L = (L Pt q) p ' q . It is easy to see that L\ U L2 = L\ U L2, and similarly for intersections and 

o 

also for the operations L, and 



n P , q ( L ) = 




We now turn to the proof. Assume first that such a shift A exists, and let L C S be a 
lower set. Then the points of L n Z p,q will be obtained by moving some of the points of 
A, and those points must come from An L. Hence A defines a bijection between a subset 
of An L and L n Z p < q , and this proves the desired inequality. We now keep p and q fixed, 
we denote by V the set of pairs (A,S) that satisfy (i) and have \S\ = n\A\, and we prove 
that (ii) holds for any (A, S) £ V by induction on The starting point is \S\ = n. 

Then A is forced to be {(0, 0)} (apply the condition to L = {(0, 0)}), S is forced to be the 
rectangle R(p,q) (apply the condition to L = S), hence A = S n Z p,q = {(0,0)}. Assume 
now that (A,S) does satisfy (i), and the implication has been shown for all pairs (A',S') 
with \S'\ < \S\. If one of the points of A can be moved one step upwards, or one step 
to the right, so that the condition (i) is still satisfied (by S and the new set A), then we 
perform the move. We repeat this if still possible. This process will stop at some point 
(e.g. because A cannot be moved outside S). Hence, we may assume that A is "maximal", 
in the sense that none of its elements can be moved any further without violating (i). We 
will show that A = S n Z p ' q . We first prove the following: 

Claim 1: If L C S is a lower subset with \A H L\ = n Ptq (L), then either L = S, or 

A n L = z p < q n L. 

Proof of the claim: We first show that (AnL,L) € V. First of all, |L| = n p>q (L) = \AC\L\. 
Secondly, for any lower set P C L one has P C L, hence LnP = LnP = P. This, and 
the fact that (A, S) does satisfy (i), imply that 

n P , q (P) = n Ptq {L n P) < \A n (L n P)\ = \(A n L) n P\ 

for all P C L lower. Hence (A n L, L) G V. Assume now that L ^ S. Then, by the 
induction hypothesis, there is a shift that moves A n L to Z p,q n L = Z p,q n L. Since this 
shift stays inside L, it does not touches the elements in A \ L, hence it can be viewed as 
a shift A of the entire A, which leaves A \ L intact. In other words, the image A' of A is 
given by 



A' n L = Z p ' q n L, A'\L = A\L. 
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Now, for any P C S lower, one has 

\A' n p\ = \a! n (P n L)\ + \A' n (P \ L)| 
= \A' n (P n L)\ + \An (P \ L)\ 
= \{A' n L) n {P n L)\ + \An (P u L)\ - \An L\ 

> n M {P n L) + n Piq (P U L) - n pq (L) 
= n P!q (R) 

In the inequality above we used that (A' n L, L) and (A, S) satisfy (i), and that \A D L\ = 
n P) q(L). This shows that (A',S) G V. From the maximality of A, the shift must be the 
identity shift, which precisely means that Af) L = Z p,q n L. This concludes the proof of 
the claim. 

Claim 2: If e G A \ (?) is an element with the property that e r ^ A, then 

e G d y {S)Ud e (S). 

Proof of the claim: Denote by A r e the new set obtained from A by moving e to e r . By 
the maximality of A, [A r e ,S) does not satisfy (i), hence we find a lower set L such that 
\A r e n L\ < n Ptq (L). But |A£ n L\ is either |An L\ - 1 (if e G L, e r L), or |A n L\ (in the 
remaining cases). Since \A (~l L| > n p ^(L), we must have 

n L| = n P) q(L), e G L, e r L. 

We now use the previous claim. Since e G (^4 n L) \ Z(p, q), we must have L = S. In turn, 

o o o 

this implies that L=S, hence, since e r ^ L, we have e r ^S- On the other hand, applying 
(i) to the largest lower set which does not contain e, we immediately see that S contains 

o o 

at least one element in Z p,q larger then e. In other words, e GS 1 . But the elements e GS 

o o o 

with e r are exactly those those situated on d y (S) U d e (S)- This concludes the proof of 
the claim. 

Clearly, one can replace e r by e u , and arrive to a similar conclusion. In particular, if 

o 

e G A \ Z p ' q , then either e r or e u must be in A (otherwise we must have e G d e (S), which 

o 

is a contradiction because the extremal boundary points of S belong to Z p ' q ). 

We are now ready to prove that A = L n Z p ' q . If this is not so, then we pick up an 
extremal element e of A \ Z p ' q . That means that e G A \ Z p ' q , and A \ Z p ' q contains no 
other element whose coordinates are greater or equal to the coordinates of e. Then at least 
one of the positions e u or e r are not in A. Otherwise, since e is extremal, e r and e u would 
both be in Z p ' q , and that is clearly impossible. Hence, combined with the previous remark, 
one (and only one) of these two elements are in A. Hence we may assume that 

e r iA, e u £A. 

From Claim 2, and the extremality of e, we deduce that 

e r ed y Cs)ud e (°s), e u eZ p > q . 

o 

Let e' be the exterior boundary point of S with the property that the segment / = [e, e'\ is 
parallel to OY (see Fig. [SJ). We apply condition (i) to the lower set S\Q e , where Q e is as 
in the picture, (i.e. S \ Q e is the largest lower subset of S which does not contain e). We 
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immediately get \l n Z p,q \ > \l fl A\. Since e is the only point of A situated on I but not on 
Z p,q , we find a point / = (a, (5) 6 A n / n i^' 9 with the property that / = (a + p, /3) is in 
I n Z p,q but it is not in A. We then consider the shift that moves / to /, and leaves the 
rest of A unchanged. The maximality of A implies again the existence of a lower set L with 
the property that \A Pi L\ = n Ptq (L), f G L, f £ L. We can now use Claim 1 above. Since 
f £ L and / £ Z p ' q , one cannot have L = S. Hence we must have A n L = Z v,q n L. On 
the other hand, since L is lower and e sits below /, it follows that e £ An L =G L, hence 
e € i? 5 ' 9 . This contradicts the choice of e, and concludes the proof of the Theorem. □ 
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Figure 5. 

Let us point out another geometric property that is suggested by examples. This time, 
we will use inverse shifts, which are defined exactly like shifts, but moving downwards and 
to the left. 

Conjecture 3. If (A, S) is almost regular with respect to (p,q) -rectangular sets of nodes, 
then there is an inverse shift of A in S, which moves A to a lower set R, and S = R p,q . 

As in Theorem 14.11 (and proven using the same ideas), this geometrical condition is 
equivalent to an algebraic condition, which says that \A Pi L\ < \R Pi L\ for all lower sets 
L. Note also that one should only allow certain type of shifts which do not violate certain 
"regularity conditions" . What "regularity condition" exactly means is still to be discovered, 
but it certainly excludes moving new elements to the axes (cf. the results of the previous 
section). Also, it is tempting to combine the last conjecture and Theorem l4. ll into a stronger 
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(unifying conjecture) which states the existence of (certain) shifts of R into SnZ p ' q , which 
move (i, j) to ({p + T)i, (q + l)j), and which, at some intermediate step, cover A. 

Note that in almost all the cases we have considered so far (e.g. no mixed derivatives, 
or lower sets of derivatives), moving A backwards to a lower set was possible in only one 
way. Hence Conjecture |3] would explain the uniqueness of S (proven by us in each case 
separately. Also, one can use this conjecture as a guide for constructing interesting examples 
(e.g. where S is not unique, see Example I5.1U|) . 

5. Examples 

In this section we present several examples that illustrate the results of the previous 
sections. For simplicity, we restrict most of the examples to the case p = q = 1, i.e. the 
case of (1, 1) -rectangular sets of nodes. One of the simplifications comes from the fact that, 
in this case, the notions of regular and almost regular coincide (cf. Theorem 13. llj) . Hence, 
unless otherwise specified, the term "regular" in this section stands for "regular with respect 
to (1, 1) -rectangular sets of nodes". Passing to the general case of (p, q) -rectangular sets of 
nodes require some care on almost regularity versus regularity, but many of the arguments 
remain the same. 

Example 5.1. Let A be as in Fig. El Then there is no lower set S which makes (A,S) 
into a regular scheme. Assume there is one. First of all, it must be i?(5,5), as implied 
by Proposition 13.11 In this case however, the inequality ® (insured by Theorem 13.11)1 is 
violated by the lower set L drawn in the picture. Hence, there is no lower set S which 
makes (A, S) regular. 

s(Z 1>1 ,A) = 0. 

Note also that after moving any of the first two points of A situated on the line y = 3 one 




x 



S\L 



points in A 



Q : points in S with 
even coordinates 



Figure 6. 



step downwards, the condition is no longer violated, and one can actually show that the 
resulting schemes are regular. However, moving the last point on that line (i.e. (3,3)) one 
step downwards, produces a scheme which is still non-regular since the same condition is 
violated (this time by L' shown in Fig. Q). 

Example 5.2. Consider now A as in Fig. |HJ As above, regularity forces S = R(5,5). 
This time however, the inequality (j3J) is satisfied. Let us show that (A, S) is regular. We 
first remark that S y (A) is (3, 3)-rectangular, hence we can use Theorem 13.21 to reduce the 
regularity of (A, S) to the regularity of several univariate schemes. In turn, the univariate 
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Figure 7. 



schemes are being taking care of by the Polya condition (see 12. 7|) . When p, q > 1 the 
same arguments apply to conclude that (A, S) is almost regular (with respect to (p, q)- 
rectangular sets of nodes) if and only if S = R(2p + 3,2q + 3). Moreover, given the set 
of nodes Z, Theorem 13.21 rephrases the regularity of (Z, A, S) in terms of the regularity of 
certain induced univariate schemes which are easier to handle. 




X : points in A 

Q : points in S with even coordinates 
■ :S 



Figure 8. 



Example 5.3. We consider now the variation of the previous example shown in Fig. |UJ 
One cannot apply Theorem 13. 21 directly, but one can first invoke Theorem l3.6l to re-arrange 
the points of A on OX to occupy the first three positions. Then S y (A) becomes rectangular, 
and Theorem l3.6l can be used. Alternatively, one can use Theorem 13.61 twice (once on each 
of the axes) to and then reduce the problem to the one of the previous example. 

Example 5.4. Let A be as in Fig. 1101 Using Theorem 13. 5[ we see that S must contain 
So shown in the picture, and it must be contained in i?(5, 5). Using Theorem 13. 1 II (namely 
that S is of type R 1 ' 1 for some lower set R), we see that S is obtained from So together with 
a copy of the rectangle R(l, 1). But there are only two ways one can add such a rectangle 
to So to obtain a lower set, and the two possibilities are shown as Case 1 and Case 2 in the 
picture. In the first case, the inequality (j3J) (insured by Theorem 13. lift with L = S\{(3, 1)} 
is not satisfied. The situation is different in the second case, when we obtain a regular 
scheme. To see this, one first uses Theorem 13.61 to remove the last point of A situated on 
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Figure 9. 



OY, and then one treats the remaining scheme as in Example 15.21 Hence, again, there is 
only one S which makes the scheme regular. 



Example 5.5. Let us now describe two examples which show that the inequalities (0} do 
not imply regularity with respect to (1, 1) -rectangular sets of nodes. Consider the scheme 
(A, S) appearing on the left hand side of Fig. 111! It does satisfy the desired condition, but 
it is not regular. To see this, we remark that A can be obtained from A (see right hand 
side of Fig. II If) by removing the last element from OX. Hence we can apply Theorem 13. 61 
to (A, S) to conclude that it is regular if [A, S) is. But this cannot happen because the 
regularity of (A, S) and the fact that A is lower would imply that S = A 1,1 (cf. Theorem 
13.41) . which is not the case. We should say here that what causes the non-regularity in this 
example is another simple condition that must be satisfied by all regular schemes (and, in 
this example, it is not): the number of points of A on the line y = 1 cannot exceed those 
on the line y = 0. 

Example 5.6. A bit more subtle is the example drawn in Fig. 1121 which still satisfies the 
inequalities Q. We advise the reader to try to guess a "general regularity condition" that 
is broken in this example. We now give an argument that proves that (A, S) cannot be 
regular. Assume it is. We first move the last two elements of A on OX to new positions to 
get the scheme (A',S') of Fig. 1121 By Theorem 13.61 (and the univariate Polya conditions, 




Case 1 (non-regular) 



Case 2 (regular) 



Figure 10. 
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by removing the last element of A on OX 

Figure 11. 

see 12 .7(1 . (A', S') is still regular. Now, choose L as in the picture. One has \L\ = 4\A n L\, 
i.e. we are in the limit case of the Polya condition. Hence (see 12.7(1 (A n L, L) must be 
regular too. But, by the same arguments as in the example above, (A n L, L) cannot be 
regular. 
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Figure 12. 

Example 5.7. Let us return to the set A appearing in the last example in 15.51 (Fig. 112(1 . 
and look for all 5's which make (A,S) regular. As in Example 15.41 there are two possible 
cases. One of them is precisely the one treated in Example 15. 51 while the other one is shown 
in Fig. 11^1 The last one is regular. To see this, one first uses Theorem IH.fil to remove the 
last point of A situated on OY. The resulting scheme {A\ S') has S X (A) rectangular, hence 
we can use Theorem 1^.21 (or, more precisely, the version obtained by interchanging x and 
y). The subsequent univariate schemes are regular (by 12.71 again). 

Example 5.8. A similar example is obtained by considering A as in Fig. 1141 As before, 
S must be obtained from So in the picture by adding one copy of R(l, 1). There are two 
possible ways to do that, but only the one shown in the picture produces a regular scheme. 
However, we do not know how to use the general results of the previous sections to prove 
the non-regularity of the other scheme. 

Example 5.9. Consider now A as in Fig. 1151 As before, Theorem 13.51 tells us what S 
must be around the axes. We then have to fit three more (disjoint) copies of R(l, 1) inside 
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Figure 14. 



R(5, 5) to get the lower set S. This time, this is possible in only one way (as in the picture). 
One sees that the inequalities (j3J) are satisfied, and one can actually show that the scheme 
is regular. 

Example 5.10. In all the previous examples, given A, there was at most one S making 
(A, S) regular. Here is an example where two such 5"s can be chosen. Consider A as in 
Fig. 1161 As in the previous examples, S is contained in R(5, 5) and must be obtained by 
adding one copy of R(l, 1) to the blow up (I2) ' of the triangle T^. This is possible in two 
ways, with the resulting S"s: S y (2, 1, 1) 1,1 and S y (2, 2, l) 1,1 (see also Fig. IT7jl . That both 
resulting schemes are regular follows again by removing one point and then using Theorem 

It is interesting to discuss this example in the light of Conjecture 01 Note that, in all 
the examples we have considered so far, moving A backwards to a lower set was possible in 
only one way, hence the conjecture suggest the uniqueness of S (proven by us in each case 
separately). In the example under discussion, the process of "moving A backwards" to a 
lower set is not unique; there are clearly (only) two ways of doing so, as shown in Fig. 1171 
the extremal element of A can be moved one step down, or one step to the left. And this 
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Figure 15. 



A: 



x — x 



X----X 



X X- 





Two different regular schemes with the same set of derivatives A 



Figure 16. 



is how we actually constructed this example. Similarly, one can find examples where the 
number of choices for S equals a given number. 
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